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Abstract 

43 . 
O !■ We present a warped extra dimension model with the custodial symmetry SU(2)l x 

qJ SU(2)fj x t/(l)x x -Pl_r based on the flavor symmetry S4 x Z2 x Z2, and the neutrinos 

are taken to be Dirac particles. At leading order, the democratic lepton mixing is 
derived exactly, and the high dimensional operators introduce corrections of order 
A c to all the three lepton mixing angles such that agreement with the experimental 
data can be achieved. The neutrino mass spectrum is predicted to be of the in- 
verted hierarchy and the second octant of 623 is preferred. We suggest the modified 



\& , democratic mixing, which is obtained by permuting the second and the third rows of 



the democratic mixing matrix, should be a good first order approximation to under- 
standing sizable #13 and the first octant of 623- The constraints on the model from 
the electroweak precision measurements are discussed. Furthermore, we investigate 
the lepton mixing patterns for all the possible residual symmetries G u and Gi in the 
neutrino and charged lepton sectors, respectively. 
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1 Introduction 

Warped extra dimensions, also known as Randall-Sundrum (RS) models, were first pro- 
posed to solve the gauge hierarchy problem in the standard model (SM) pQ. In addition, 
they provide a novel and powerful framework for understanding flavor physics. The ob- 
served SM charged fermion mass hierarchies can be naturally generated due to the overlap 
of the fermion and Higgs wavefunctions, if the SM fermions and gauge bosons are allowed 
to propagate in the bulk [2]. Indeed, from this approach, both the quark mass hierarchies 
and the CKM mixing angles can be accurately reproduced as shown in Refs. [2111], and 
the Yukawa couplings can be naturally of order one with a completely random pattern. 
However, the electroweak precision data mainly from the Peskin-Takeuchi T parameter 
turn out to be so restrictive that the masses of the Kaluza-Klein (KK) modes of bulk 
fields are pushed to somewhat high scales that make their discovery at the LHC extremely 
challenging [5] . To reduce the size of the oblique corrections, one can extend the SM gauge 
symmetry to a bulk custodial symmetry [HIE], introduce large brane kinetic terms [5] or 
promote the Higgs to a bulk field [9]. 

In general, the 5D fermion bulk mass parameters, which determine the localization of 
the zero mode field in the extra dimension, are not degenerate for different fermions in 
order to account for the mass hierarchy. As a result, dangerously large flavor changing 
neutral current (FCNC) processes arise already at the tree level through the KK gauge 
boson exchange [T0HT6] . To soften the FCNC constraints, one could introduce additional 
bulk flavor symmetry such as the 5D minimal flavor violation [T7J[T8] or some dynamical 
mechanisms [19J. Another interesting way is to impose a discrete flavor symmetry to 
provide extra flavor protection [20], no additional gauge bosons have to be introduced in 
this setup, and we could possibly obtain realistic lepton masses and flavor mixing pattern 
through the spontaneous breaking of the discrete flavor symmetry. The idea of combining 
the discrete flavor symmetry and extra dimension is quite attractive and has already been 
discussed in the literature within the context of large extra dimensions [21} [22], warped 
extra dimensions [23] and holographic composite Higgs models [21]. However, all these 
models try to generate tri-bimaximal neutrino mixing [25] at leading order (LO), where 
#i3 = 0°, #i2 — 35.3° and #23 = 45°. It is possible to generate non-zero #13 by the next to 
leading order (NLO) corrections, but in the generic case, the induced deviations should 
be of the same order for all the three mixing angles. Due to the strong constraint from 
the accurately measured #12, these models are expected to give rise to #13 at most of order 
A^, where A c ~ 0.23 is the Cabibbo angle. Recently the neutrino oscillation experiments 
T2K [26], MINOS [21], Double CHOOZ [2H], Daya Bay [29] and RENO [30], together with 
the global fitting of mixing parameters [3TH33J . confirm a sizeable reactor angle #13 ~ 9°. 
Consequently, the tri-bimaximal mixing as a LO approximation is strongly disfavored, 
unless the underlying theory is capable of providing sufficiently large corrections to 9 Yi 
without affecting too much the solar angle. 

Flavor symmetry has been widely applied to explaining the structure of the leptonic 
mixing angles in the past. Roughly speaking, there are two approaches based on flavor 
symmetry to understand the current lepton mixing parameters, particularly the sizable 
6*i3 and non- maximal #23- The first one is to construct models that give some new mixing 



textures, which are in good agreement with the present data and initially admit a non-zero 
9 13 . The so-called Toorop-Feruglio-Hagedorn (TFH) mixing pattern, which gives sin 2 6* 13 = 
(2 - v/3)/6 ~ 0.045, sin 2 12 = (8 - 2 v / 3)/13 ~ 0.349 and sin 2 9 23 = (5 + 2^)/ 13 ~ 0.651, 
is a viable substitute for the tri-bimaximal mixing, and it can be naturally produced from 
the A(96) family symmetry group J31H3S]- Under the hypothesis of Majorana neutrinos [J, 
an extensive general scan of all finite discrete groups with order less than 1536 is performed 
to obtain the corresponding predictions for lepton mixing angles [37]. Only three groups 
A(600), (Z 18 x Zq) xi S3 and A(1536) are found to predict lepton mixing angles within 
3er range of the current global fits. That is to say, the order of the flavor symmetry 
group is required to be rather large □ and the resulting flavor models become typically 
complex [32] • The second approach is starting from some texture which is outside the 
current 3<r range, but the sizable NLO corrections could pull the three mixing angles into 
the experimentally fovared range. In Ref. [40], the well-known bimaximal mixing [41] 
with #13 = 0, 9 12 = 6*23 = 45° was taken as the first order approximation, and the model 
was tactfully constructed such that # 12 and 9 13 were corrected by terms of order A c while 
#23 was unchanged at this order. As a result, the mixing angles could be in accordance 
with the measured values. In this work, we shall investigate the democratic (DC) mixing 
pattern [42] as the LO mixing pattern. In order to be compatible with experimental 
data, all the three mixing angles should undergo corrections of order 0.1 ~ 0.2, i.e. the 
NLO corrections should be of order A c , as has been pointed out in Refs. [43][44]. The 
corresponding models would be somewhat more natural and more easily constructed that 
the bimaximal case, since we don't need specific dynamical tricks to realize the different 
order of magnitudes of the NLO corrections to the three mixing angles. 

It is now generally accepted that neutrinos have tiny masses. However, the nature 
(Dirac vs. Majorana) of neutrinos and thus the origin of neutrino masses remain unknown. 
Historically in the 4-dimensional (4D) space-time, the preference is that they are of the 
Majorana type, since the smallness of the neutrino masses can be elegantly explained by 
the see-saw mechanism [45]. The see-saw mechanism is also valid in the 5-dimensional 
(5D) case [45JII7]. Moreover, warped extra dimension can generate small Dirac neutrino 
masses without fine-tuning from the appropriate localization of neutrinos in the bulk [45] . 
In addition, it was argued in [49], that in RS type scenarios, constraints from lepton flavour 
violation favor Dirac neutrinos over Majorana neutrinos. Therefore we restrict ourselves 
to the case of Dirac neutrinos. Another motivation for considering Dirac neutrinos is 
below: following the bottom-up method presented in Ref. [5THI5TJ , we find that there is no 
discrete group that contains all of the symmetries needed for the DC mixing, if neutrinos 
are Majorana particles. Within the framework of warped extra dimensions, we shall show 
that the DC mixing pattern can be naturally produced at LO from £4 flavor symmetry. 



lr They assume that the remnant symmetries in the neutrino and charged lepton sectors are Zi X Z2 
and Z3, respectively, and the left-handed lepton doublets are assigned to the three dimensional irreducible 
representation of the flavor symmetry group. 

2 If only Z2 instead of Z2 X Z2 subgroup is preserved in the neutrino sector, the leptonic mixing angles 
can not be completely fixed, although they are usually correlated with each other. One can tune the 
parameters, which are undetermined by the flavor symmetry, to account for the data on neutrino mixing. 
In this case, the order of the flavor symmetry group can be smaller [38] . 



The model is built in such a way that the DC mixing receives corrections of order A c from 
the higher order terms. As a result, all the three leptonic mixing angles are corrected by 
terms of order A c , and thus, an agreement with the experimental data is achieved. 

The paper is organized as follows. In section [2] we investigate the general properties of 
the Dirac neutrino mass matrix leading to the DC mixing in the flavor basis. In section 
[3] we construct our model in the context of custodial warped extra dimensions based on 
£4 flavor symmetry, where the Higgs field is assumed to propagate in the bulk and the 
DC mixing pattern is realized in the LO. In section H] the corrections to the DC mixing 
angles from higher dimensional operators are analyzed. In section [5] the constraints from 
electroweak precision measurements are studied. In section [6] we explore all the possible 
leptonic mixing patterns within S4 by scanning the possible G u and Gi subgroups preserved 
in the neutrino and charged lepton sectors, respectively. Finally section [7] is devoted to 
our conclusion and summary. 



2 Democratic mixing 

In a specific phase convention, the democratic (DC) mixing matrix is given by 
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It gives rise to the three leptonic mixing angles as follows 
0° c = 45°, 0£° = arctan V2 ~ 54.7°, 
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The DC mixing is closely related to the well-known tri-bimaximal mixing [25] , the former 
can be obtained by transposing the latter one and then permuting its rows and columns 
respectively: 
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where Utb is the tri-bimaximal mixing matrix 
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In the basis where the charged lepton mass matrix is diagonal, the hermitian combination 
of the neutrino mass matrix M. v = m v v(V v for the DC mixing is given by 
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where 777.1,2,3 are the light neutrino masses. As a result, the most general mass matrix M. v 
corresponding to DC mixing is of the form 

(x iy —i\/ r 2y \ 

-iy z (z-x)/V2 , (6) 

iV2y (z-x)/V2 (x + z)/2 J 

where the parameters x, y and z are real. The squared light neutrino masses can be 
expressed in terms of x, y and z as m\ = x — v3 y, m\ = x + a/3 y and m\ = (3z — x)/2. If 
we perform a unitary transformation i>l — > Gvl, where ul = (u e , u^, Vtjl is the left-handed 
neutrino field, then hA v transforms as M. v —} G^M. U G. The hermitian combination M. v 
in Eq. (jHj) or Eq. ([6]) for the DC mixing would be invariant, if the transformation G takes 
the form 
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where a and /3 are arbitrary real parameters. It is straightforward to check that the above 
symmetry transformations satisfy the following multiplication rules 

G x (ai, 0i)Gi(a 2 , /3 2 ) = Gi(a x + a 2 , fix + /3 2 ) , 
G 2 (ax, f3 1 )G 2 (a 2 , (3 2 ) = G 2 (ai + a 2 , f3 x + (3 2 ) , 
G l (a l ,f3 1 )G 2 (a 2 ,p 2 ) = G 2 (a 2 ,(3 2 )G l (a u (3 1 ). (9) 

Therefore, the symmetry group of Ai u , which is generated by the symmetry transforma- 
tions Gi(a,(3) and G 2 (a,/3), is isomorphic to 17(1) x U(l) x f/(l) x 17(1) = Q v . Generally 
the symmetry group of the Dirac neutrino mass matrix for mass independent mixing pat- 
terns is Q v = U(l) x 17(1) x 17(1) x U(l) in the charged lepton diagonal base, although 
the specific forms of the symmetry transformations Gi(a,(3) and G 2 (a,(3) depend on the 
mixing texture. On the other hand, if the hermitian combination of the neutrino mass 
matrix M. v = m v m v is invariant under the Q v symmetry group generated by Gi(a,/3) 
and G 2 (a>,(3) in Eqs. (J7J) and (jSJ), it would be diagonalized by the DC mixing matrix 
precisely. In fact, it is sufficient to require M, v invariant under one element h G Q v 
with non-degenerate eigenvalues for M. u being diagonalized by DC mixing. For example, 
M. u = m u m\, invariant under an order three element h with 

/i = G 1 (47r/3,0)G 2 (2vr/3,0) = i I 1 1 \/2 | , (10) 

would be diagonalized by the DC mixing matrix. We note that this unitary transformation 
of h is exactly the representation matrix for the S4 generator T in three dimensional 
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Table 1: Representation assignments for the bulk fields and the localized scalars (i.e. flavon fields) under 
the gauge symmetry SU(2)l x SU(2)r x U(1)x x Plr and the flavor symmetry group Si x Z 2 x Z' 2 , 
where the superscripts are the family indices. 



irreducible representations, as presented in the Appendix. This is the reason why we can 
reproduce DC mixing exactly at LO from the finite £4 family symmetry group, as will be 
shown in the following. It is not necessary that the flavor symmetry group is large enough 
to include the left-handed (LH) neutrino symmetry group Q v = U(l) x £7(1) x U(l) x U(l) 
as a subgroup so that it can be preserved in the neutrino sector. On the contrary, if 
neutrinos are Majorana particles, generally the neutrino mass matrix symmetry group 
K4 = Z 2 x Z 2 should be a subgroup of the full flavor symmetry group, or one Z 2 belongs 
to the flavor symmetry group while the other is accidental. 

3 The structure of the model 

In this section, we present a warped extra dimension model of leptons based on 5*4 x 
Z 2 x Z 2 flavor symmetry. The LO lepton mixing is of the DC form in this model, and its 
phenomenological implications will be discussed. 



3.1 The model 

Our model is formulated on a slice of AdSs space. The background metric is of the 
usual Randall-Sundrum form [1] 



ds 2 = e-^^dx^dx" - dy 7 
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where r]^ = diag(l, —1, —1, —1), the fifth coordinate is restricted to the interval < y < L, 
and k ~ 0(Mp{) is the AdSs curvature scale. As usual the UV (Planck) brane is localized 
at y — 0, and the IR (TeV) brane is at y = L. It is useful to define a parameter k' = ke~ kL 
for the following discussions. To solve the gauge hierarchy problem, k' has to be at the 



TeV scale. The SM electroweak gauge group is extended to 

G bulk = SU(2) L x SU(2) R x U(l) x x P LR , (12) 

in the bulk. The 5D gauge coupling constants of SU(2) L , SU{2) R and U(l)x are g 5 L, g 5R 
and g^Xi respectively. The presence of SU(2) R implies an unbroken custodial symmetry 
which eliminates the excessively large contributions to the T parameter due to the Kaluza- 
Klein (KK) excitations of the bulk fields. The discrete symmetry Plr, interchanging the 
two SU(2) groups, is proposed to protect the SM Z boson couplings to the LH down-type 
quarks, in particular the ZhiJ)^ coupling, from large corrections [52]. Clearly the P^ R 
symmetry implies that the gauge couplings of SU(2)l and SU(2) R have to be equal g$i = 
g 5R = g 5 . The bulk gauge symmetry Gbuik is broken to the SM gauge group SU(2)l x U{1)y 
on the UV brane by means of the boundary conditions of the electroweak (EW) gauge 



bosons. The resulting 5D U{1)y gauge coupling constant is given by g' 5 = g 5X / 'yyf + g\ x . 
In order to realize the standard EW symmetry breaking, we introduce a bulk Higgs field 
which is a self-dual bidoublet of SU{2) L x SU{2) R and can be parameterized as: 



1 / v / 2tt + -h° + iir 



H= 2{h° + i* V2n- ) ■ (13) 

It is a singlet under U(l)x with Q X (H) = 0. In the notation of Eq. (fT3|) . SU{2) L acts 
vertically while SU(2) R acts horizontally, 

H^U L HIJT, U L eSU(2) L , U R eSU(2) R . (14) 

The Goldstone bosons eaten by the gauge bosons are denoted by 7r =l= and tt°. h° is the 
physical Higgs boson, and it develops a 4D effective vacuum expectation value (VEV) near 
the IR brane, which eventually leads to EW symmetry breaking. In order to yield a light 
Higgs zero mode, the whole bidoublet H has to obey the (++) boundary conditions (BCs), 
where "+" stands for a Neumann boundary condition. Performing the KK decomposition 
of the Higgs field H(x' J ', y), we have 

H(x fl , y) = - 1 =H(x fi )h(y) + heavy KK modes , (15) 

and we can choose the BCs such that its profile h(y) takes the simple form [55] lf| 



S§«'-«--'). (16) 



The parameter (3 is related to the Higgs bulk mass m^ by (3 = y/A + mjj/k 2 , which 
controls the localization of the profile of the VEV in the bulk. In order to localize the 
Higgs VEV close to the IR brane such that the electroweak symmetry is broken at the TeV 



3 The bulk equations of motion admit two solutions for the profile of the bulk Higgs: one is proportional 
to e( 2 ~ /3 ) fc ( y-L ) and the other is proportional to e ( 2+/3 ) fc ^~ L ). We could choose the BCs to select one of 
the above two solutions. In the present work, we take the first one in order to generate light neutrino 
masses of correct order of magnitude, as will be discussed in the following. 



scale, the bulk mass squared m 2 H should be negative but it is bounded by m 2 H > —4k 2 from 
both the Breitenlohner-Freedman bound [53] and the requirement that the dual theory 
resolves the gauge hierarchy problem [51] . Furthermore, the VEV of the Higgs zero mode 
is of the form 

«=(,;f). <-> 

where v denotes the effective 4D VEV of the zero mode of h°. As has been shown above, a 
novel feature of the bulk Higgs is that the Higgs VEV has a profile which extends into the 
bulk and is no longer a constant. This provides new model-building possibilities, and it is 
relevant to neutrino mass generation in the present work. Moreover, a number of existing 
studies on bulk Higgs scenarios [55j|56] demonstrate that the bulk Higgs can improve the 
naturalness of the warped models and lead to better agreement with the experiment data. 
The matter fields of the model transform in the non-minimal representations under 
SU{2) L x SU{2) R [7J|57]. There are three multiplets per generation {i = 1,2,3): 

(£l) OQ = { X e, { _ + )Q *«(++)_ J ~ (2 ' 2) °' ^ - V ( -" )0 ~ (1 ' 1)0 ' 
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where SU(2) L indices are denoted by Latin letters, and SU(2) R indices are denoted by 
Greek letters. The subscripts of the various components indicate their electric charges, 
and the U(\)x charges of all the above leptonic matter fields are vanishing. Note that in 
the notation of Eq. (fT8j) . the matter fields transform in the following way under the gauge 
transformations Ul E SU{2)l and U R G SU(2) r , 

&^U L &Ul, Ti^U L Tpl, T A ^U R TIU R . (19) 

Here we have denoted Tj = r a (T 3 4 ) a and T\ = T a (T 4 4 ) a , where r a = a a /2 and r a = a a /2 
are the generators of the fundamental SU{2) L and SU{2) R representations respectively 
with a a and a a being the Pauli matrices. The signs in brackets in Eq. (fT8l) indicate the 
BCs for the LH fermion mode on the UV and IR branes, respectively, where " + " denotes 
a modified Neumann BC and " — " stands for a Dirichlet BC. The corresponding right- 
handed (RH) modes, which are necessarily present in a 5D theory, obey opposite BCs. As 
only fields with (++) BCs have massless zero modes, up to small mixing effects with other 
massive modes due to the transformation to mass eigenstates, the low energy spectrum 
will contain a LH doublet {tl , t^ ) and two RH singlets v R and L^ ' (i = 1, 2, 3) for 
each lepton generation. The field contents of the SM are reproduced exactly. 

The full flavor symmetry of the model is 5*4 x Z2 x Z' 2 . The three generations £\ (i = 
1,2,3), which contain LH zero modes (£ u £ Al ) transforming as an SU(2) L doublets, 
are embedded into an S4 triplet 3. The three generations of ££ (^ = 1)2,3) containing 
RH neutral singlets u R , are also assigned to form an S4 triplet 3. While each generation 
of £3 (i = 1,2,3), which allows a RH singlet L^ with electric charge —1, transforms as 



nonequivalent one-dimensional representation 1 or V Q The auxiliary symmetry Z 2 is 
imposed to further distinguish the three generation £ 3 ' ' fields, and the purpose of the 
other Z 2 is to eliminate certain dangerous operators that would otherwise contribute to 
the neutrino masses at LO. 

The flavor symmetry is broken by a set of flavon fields, which are all singlets under 
the bulk gauge group and localized on the two branes. Concretely, the flavons \ and 
ip, transforming as S4 triplets 3 and 3' on the IR brane, respectively, are responsible for 
the flavor symmetry breaking in the charged lepton sector. The scalars ( ~ 1, ~ 3' 
and p ~ 2 are localized on the UV brane, and they break the S4 flavor symmetry in 
the neutrino sector. The field contents of the model and their transformation properties 
under the gauge and flavor symmetry groups are summarized in Table [TJ By putting the 
64 breaking fields \ an d <£> on the IR brane, and the remaining flavons £, <fi and p on the 
UV brane, the questions regarding vacuum alignment are eliminated, as the two sets of 
flavons are completely sequestered [21]. As usually done in such models, we shall assume 
that the flavor symmetry S4 x Z 2 x Z' 2 is spontaneously broken by the VEVs of the flavons 
on the branes as follows 

(x) = (o,o,ik, (</?> = (0,1, OK, 

(0 = v o (</>) = (0,^,1^, (p) = (l,0)v (20) 

In the chosen basis which is presented in the Appendix, it is straightforward to check that 
we have TST( X ) = S 2 ( X ) = TSTS 2 ( X ) = ( X ) and TST(<p} = -S 2 (<p) = -TSTS 2 (<p) = 
(if). Therefore the VEV of \ breaks the flavor symmetry 5*4 to its subgroup K4 = 
{1, TST, S 2 , TSTS 2 }, while the VEV of (p breaks S4 to Z 2 . On the other hand, these are 
the most general vacuum configurations which preserve the above K4 and Z 2 subgroups, 
respectively. In the same manner, we have T(4>) = (</>) and likewise it is the most general 
VEVs invariant under the action of T. Consequently, the vacuum alignment of ( and 
breaks S4 to its Z 3 subgroup generated by T. 

In the warped extra dimension, the fermion mass hierarchies are naturally explained 
through wavefunction localization, which is determined by the bulk fermion mass. The 
bulk mass Lagrangian for the matter fields ££, £,\ and £3, invariant under the flavor sym- 
metry, reads 

I- 3 3 

+ Cfl k [(Tl) a (T 2 ) a + (Tl) a (T 2 ) a ] + c T k [rt) a (T|) a + (TlUTi) a ] } , (21) 

where G = e~ 8ky is the determinant of the RS-metric, the summation over repeated indices 
is understood, and again SU(2) L indices are denoted by Latin letters while SU(2) R indices 
are denoted by Greek letters. The bulk Dirac masses are conveniently parameterized in 
terms of the AdS curvature k. In general, the bulk mass would be a hermitian 3x3 matrix. 



A Si group has only two nonequivalent one-dimensional representations 1 and 1', consequently S4 can 
not differentiate the three fields £ 3 : ' completely. 



Once we impose the S4 flavor symmetry in the bulk of the theory, the three generations of 
the fermion bidoublets £\ will share one common c— parameter which we label by q, since 
they are embedded into a triplet 3 of S4. In the same way, the 5*4 flavor symmetry implies 
a family independent bulk mass parameter c v for the three gauge singlets £\, which are 
also assigned to a triplet 3. The three generations of £3 = T\ © T\ transform as 1, 1' and 
1, respectively, and thus there are three separate c— parameters c e , c M and c T for these 
fields. Note that the Plr symmetry imposes a common bulk mass parameter for T\ and 
T4 within the same generation. 

The most general Yukawa interactions compatible with both gauge and flavor symme- 
tries, can be written as 
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where the dots stand for higher dimensional operators that will be addressed later in the 
paper, and all the above operators should be contracted into £4 singlet with the help of 
Clebsch-Gordon coefficients listed in the Appendix. All the SM fields would be identified 
with the 4D components of the zero modes in the KK decomposition of the bulk fields 
in the following, and this is the so-called zero mode approximation (ZMA). In general, 
the mixings between the zero modes of the leptons and their KK excitations could be 
induced by the Yukawa couplings in Eqs. (122J1 and (1231) . However, such a mixing effect 
is not significant, and it is at most of order m T jmxK ~ 10~ 3 in rough estimate. Many 
phenomenological analyses have demonstrated that the effects of the fermionic KK modes 
are really subleading and therefore can be neglected [52j|58] • Once the electroweak and £4 
flavor symmetries are spontaneously broken, both neutrinos and charged leptons acquire 
masses via the above Yukawa interactions. The size of the lepton masses is determined 
by the wavefunction overlap of lepton zero modes with opposite chirality and the VEV 
profile of the Higgs field. The IR boundary terms in Eq. (|22|) . which connect £} with £3, 
are responsible for the charged lepton masses. Inserting the flavon VEVs shown in Eq. 
f)20p . we find that the charged lepton mass matrix is exactly diagonal with 



m P 



m, 



m T 



Uc 



v x v v 



(ALf 2 A 2 V2 



F(y,/3,c h c e 



Vii 



Vy v 



(AL) 3/2 A v/2 
Vt v x V 



F(y,^,ci,c IJ ,] 
F(y,P,ci,Cr) 



(AL) 3/2 A v/2 

where we introduce the overlap function F(y,/3, cl, Cr) for convenience 
F(y, ?, c L , c R ) = e- k yh(y)f { °\y, c L )f R °\y, c R ) 

(fcL) 3/2 e -^-P)kL e (2-f3-c L +c R )ky 



(24) 



2 (1 - P) (1 - 2c L ) (1 + 2c R ) 



[1 _ e -2(l-/3)fcL] [ e (l-2 Ci )fci _ X] [ e (l+2 Cis )fcL _ ]_] 



(25) 
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-(0) 



Here f L (y, cl) is the profile of the LH fermionic zero mode, with the form [T0tH8 



It exists only for (++) BCs for the LH mode. f R (y,cn) is the profile of the RH zero 
mode, which is given by 



/f^W l^fE^ *- < 27 > 



Here we have performed KK decompositions for fermions as 

-ky °° 



v 7 ! 

v n=0 



Then straightforwardly, we have 



Fly = L 8 c L c R ) = (kLf\ 2 (1 - fl (1 - 2c £ ) (1 + 2c fl ) ^ CL+CR)kL 

*\y L,,p, C L , C R ) yKL,) y ^ _ e -2(l-/3)fcL] [ e (l-2c L )fcL _ ]_] [ e (l+2 CR )fcL _ 1] 

which weakly depends on the Higgs localization parameter /3. The charged lepton mass 
hierarchies are mainly governed by the exponential factor e^~ CL+CR ' kL j n the usual way. 
Note that the flavor symmetry breaking provides an additional suppression factor v^/A 
(v x /A) for the electron mass with respect to the tau (muon) case. 

The UV terms in Eq. (I23J1 connect £* with ^ an d are responsible for the neutrino 
Yukawa coupling. It is an interesting new feature due to the bulk Higgs introduced in the 
present work, implying that the Dirac neutrino masses can only come from the IR terms if 
the Higgs field is confined on the IR brane. With the vacuum in Eg. (120)) . the Dirac mass 
matrix for the zero mode neutrinos is given by 



a 



b ~^[2V 



m v =\ -b a I 4= , (29) 

.V260 a J V2 

where 

o=^-^F(i/ = 0, J 9 l c l ,(v), b = ^—^F(y = 0,^c l ,c I/ ). (30) 

In this context, the tiny neutrino masses are dominantly determined by the wavefunction 
overlaps of the Higgs and the UV scalar ( or <fi. Since the Higgs profile is localized 
close to the IR brane, the neutrino masses are highly suppressed. The suppression factor 
F(y = 0, a, C£, c v ) is given by 



F(y = 0, ft c„c„) = (kLf^ {i _ e J^-^:^]^l )tL _ T] e-"-«-.(31) 

Therefore the light neutrino masses are roughly suppressed by the factor e~^~^ kL with 
respect to the electroweak scale, and the correct scale of the neutrino masses could be 
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easily obtained by varying the Higgs localization parameter (3, as will be investigated in 
detail in the end of this section |f|. Now we come to the diagonalization of the neutrino 
mass matrix in Eq. f )29|) . which can be diagonalized by bi-unitary transformations, i.e., 

U ] L m u U R = diag(mi, m 2 , m 3 ) , (32) 



where the unitary matrices Ul and Ur are 



U L = U R 



( -i- —-i- \ 

V2 V2 U » 

i i_ fi 

Ve Ve V 3 

, JL J_ J_ i 

\ V3 v^ v^ / 



(33) 



It is exactly the DC mixing matrix. Since the charged lepton mass matrix is diagonal, 
the lepton mixing entirely comes from the neutrino sector, it is of the DC form. The light 
neutrino masses 7721,2,3 are given by a 

(\ 13 / \ 13 11 

a + iv3b) —= , m 2 — (a — iv3b) —j= , m 3 = cl—j= • (34) 

Obviously the light neutrino masses satisfy the sum rule 

m\ + m 2 = 2m 3 , (35) 

where rrii (i = 1,2,3) are generally complex numbers, so this relation cannot be used 
to exactly predict one physical neutrino mass in terms of the other two. Note that this 
sum rule is different from previous neutrino mass sum rules obtained in the context of 
Majorana neutrinos [59] . 



3.2 Phenomenological implications 

As shown above, the neutrino part of the model depends on two complex parameters a 
and b at LO. Absorbing the overall phase by field redefinition, we are left with three real 
parameters, which can be chosen to be \a\, \b\ and $ = arg(6) — arg(a) with $ being the 
relative phase between the complex numbers a and b. Now we turn to its predictions for 
the neutrino mass spectrum. From Eq. ( }34"1) we have 

,2 „,2 



mi| 2 , |m 2 | 2 = (\a\ 2 + 3\b\ 2 T 2V3|o||6| sin $) V — , |m 3 | 2 = |a| 2 y, (36) 



5 Note that if we wanted to choose the second solution for the Higgs profile which is proportional to 
e {2+0)k(y-L) ^ a vj we wou i c [ have to do is to flip the sign of j3 in all the equations. In other words, taking 
P < formally corresponds to the other solution. Then the neutrino mass scale would be suppressed by 
e -(!+^) fcL relative to the electroweak scale, generally the resulting neutrino masses would be too small to 
accommodate the measured solar and atmospheric mass-squared splittings without strong fine-tuning. 

6 As the order of the neutrino mass is less constrained experimentally so far, we can also reorder the 
light neutrino masses such that mi = (a — iv3fe) v/y/2, 1112 = (a + zv3&) v/y2 and 1TI3 = av/y/2. The 
associated unitary transformations can be obtained from Ul and Ur by permuting its first and second 
columns. The resulting lepton mixing matrix is still of the DC form, and the sign difference with respect 
to the DC mixing matrix of Eq. (TjQ) in the first row can be eliminated by redefining the lepton fields. In 
general, if the LO mixing ansatz gives rise to #12 = 45°, the lepton flavor mixing is invariant under the 
exchange of the first and the second light neutrino mass eigenstates. The bimaximal mixing is another 
typical example. 
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which leads to LI 

\m2\ 2 > (|mi 

where we have employed the constraint \m 2 \ 2 > 
experiments in the first inequality. Therefore the light neutrino mass spectrum can only 
be inverted hierarchy in the present model. Experimentally, the solar and atmospheric 
mass-squared splittings have been measured precisely. For the inverted hierarchy, they 



2 +|m 2 | 2 )/2>|m 3 | 2 , (37) 

mi\ 2 from the solar neutrino oscillation 



arc 



Am 2 



sol = K| 2 - K| 2 = (7.50±8:i|) x 10~ 5 eV 2 



Ami = M 2 - Kl 2 = (2.4271852) * lO^eV 3 • 



(3J 



We can express the three light neutrino masses in terms of r = Am 2 so J Am 2 tm , Am 2 tm and 
sin $ as 

r H * — Am 



mi 



^2 



% 



8(2-r)sin 2 $_ 
Am 



atm ' 



8(2-r)sin 2 $_ 



atm i 



Am 



atm ' 



(39) 



8(2 -r) sin 2 $ 

Note that these relations are also satisfied even if we switch the first and second light 
neutrino masses shown in Eq. ([3"4"]1 . Thus we obtain the following limit for the lightest 
neutrino mass m 3 , 

\m 3 \ >3.8 x l(T 4 eV, (40) 

where the lower bound corresponds to | sin $ | = 1 . This implies that the light neutrino 
mass spectrum has a partial hierarchy: 



mi 



4.8 x 10 _2 eV. 



\m 2 \ 



4.9 x 10 _2 eV, 



\m 3 \ 



3.8 x 10" 4 eV 



(41) 



and the sum of neutrino masses is about 9.8 x 10~ 2 eV. If sin$ is accidentally small, the 
neutrino spectrum becomes degenerate. We have plotted the neutrino masses predicted 
by the model in Fig. [U where the horizontal lines are the cosmological bounds jSU]. The 
first one at 0.60 eV corresponds to the combination of the Cosmic Microwave Background 
(CMB) anisotropy data (from WMAP 5y [61], Arcminute Cosmology Bolometer Array 
Receiver (ACBAR) [62], Very Small Array (VSA) [63], Cosmic Background Imager (CBI) 
[61] and BOOMERANG [63] experiments) plus the large-scale structure (LSS) information 
on galaxy clustering (from the Luminous Red Galaxies Sloan Digital Sky Survey (SDSS) 
[66] ) plus the Hubble Space Telescope (HST) plus the luminosity distance SN-Ia data 
of [67] and finally plus the BAO data from [68]. The second one at 0.19 eV corresponds to 
all the previous data combined to the small scale primordial spectrum from Lyman-alpha 
(Lya) forest clouds [69] • We see that the sum of light neutrino masses in the present model 
is below the cosmological bound except for extremely small sin$. 



7 In this case, one can easily deduce sin$ > from the constraint Am^ ol 
2y/3 v 2 |a||6|sin$ > 0, while one has sin$ < for the other mass ordering mi 
and 77i2 = (a + iyibj v/y2. 



\m 2 \ - |mi| = 
(a-iVSb)v/V2 
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Figure 1: The lightest neutrino mass |ttt-3 | and the sum of neutrino masses X)i=i l m »l versus | sin$|. The 
bands have been obtained by varying both &m 2 sol and Am 2 tm in their 3c experimental ranges. 



In the following, we give an illustrative example, a set of parameters reproducing the 
charged lepton masses and the neutrino mass-squared differences. We take the fundamen- 
tal 5D scale to be k ~ A ~ M Ph where M Pl ~ 2.44 x 10 18 GeV is the reduced Planck 
mass. We set the value k' = ke~ kL ~ 1.5 TeV in order to resolve the hierarchy between 
the Planck and electroweak symmetry breaking scales, and the resulting first KK gauge 
bosons with mass ttikk = 3 ~ 4 TeV are in the reach of LHC experiments. The Higgs 
VEV turns out to be v ~ 250 GeV. It is obtained by matching the precisely measured 
fine-structure constant, Fermi constant and the Z boson mass with the corresponding ones 
of the present warped model, as will be shown in section [5j The parameter V/A is fixed 
at the indicative value of 0.15 which is implied by the higher order corrections to achieve 
agreement with the measured lepton mixing angles, where V stands for a generic flavon 
VEV (v x , Dip, V£, v,/, or v p ). The Higgs localization parameter /3 will be set to 0.25 for 
demonstration in the following. 

In the charged lepton sector, we make the choice q = 0.51, c e = —0.70, c M = —0.58, 
c T = —0.502, y e = 1.248, y^ = 0.932 and y T = 2.363. These values give the charged lepton 
masses m e ~ 0.490 MeV, m^ ~ 103.3 MeV and m T ~ 1.757 GeV, which are consistent 
with the corresponding values around ITeV scale [70]. In the neutrino sector, clearly the 
light neutrino mass scale is characterized by m = F(y = 0, (3, ci, c u )Vv/(y/2A) . For the 
Higgs profile of Eq. ( 1T61) . m^ increases with the Higgs bulk mass parameter /3 sharply, 
as is displayed in Fig. |2j From this figure, we can see that the curve moves from the 
upper left corner to the lower right corner with the increase of c u . If the RH neutrino 
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Figure 2: The light neutrino mass scale mo versus the bulk Higgs parameter j3 for typical UV localized 
RH neutrino c„ = —1.2, c„ = —0.6 and IR localized RH neutrino c„ = —0.4, c„ = —0.2, where ci is taken 
to be 0.51. The horizontal lines correspond to too being 10 -4 eV and 10 _1 eV, respectively. 



is localized close to the UV brane (c u < —0.5), rriQ weakly depends on c v . Taking for 
example c v = —0.6, the neutrino mass in the range of 10~ 4 ~ 10 _1 eV can be obtained 
for /3 in the interval 0.11 < /3 < 0.31. We note that this (3 range varies slightly with 
the UV type c„. On the other hand, if the RH neutrino is localized near the IR brane 
(c v > —0.5), the localization parameter /3 should be larger to account for the correct 
magnitude of neutrino masses, since they are approximately suppressed by the factor 
e -(3/2-p+c u )kL j n ^jg case |s| yy e uge fa e value of q above, c v = —0.6, y\ = 0.95 and 

y2 = 2.49 + 0.09i. These parameters give absolute neutrino masses mi ~ 0.0496 eV, 
m 2 — 0.0503 eV and 7713 ~ 0.0107 eV. These correspond to the mass squared differences 



sol 



7.57 x 10" 5 eV 2 and Am 2 atm ~ 



2.418 x 10 3 eV 2 , which are in good agreement 



Am 

with the experimental values [31TI33] . 

8 If we had selected the other solution for the Higgs profile which is proportional to eS. 2 +P) k (v— L ) ^ rrio 
should have decreased with increasing /3, and too should have been typically of order 10 -6 ~ 10 -5 eV for 
small j3 ~ and —c v ~0(l). To produce the correct order of magnitude of the neutrino masses, the 
couplings yi and 1/2 would be very large, about several thousands, or the RH neutrinos would be localized 
extremely close to UV brane, such that their bulk mass parameters should be around tens of thousands. 
This is exactly the reason why we choose the Higgs profile of Eq. (TT6]) . 
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4 Higher order corrections 

The results of the previous section hold to the LO approximation. High dimensional 
operators, suppressed by additional powers of the cut-off A, can be added to the leading 
terms in Eqs. (122 j) and (I23p . Another motivation for considering higher order corrections 
is to obtain an acceptable lepton mixing pattern in particular to get an appreciable #13, 
since the LO exact DC mixing implies three mixing angles sin 2 9® 2 C = 1/2, sin 2 9 23 c — 2/3 
and sin 2 0^ = 0, which are sizeably different from the experimentally measured values. 

We first focus on the IR brane. The most general higher order operators invariant 
under the imposed symmetry take the following form 

*c™= E wkm l Tr ( ? ^ g > + Tr & g7 ? r > ] (*V). 

m>2,n>0 
m>0,n>l 

+ E A(m+ L 5 / 2 ) [ T r gi^) + Tr (f^Tf) ] (xV ra+1 ) 3 + h-c , (42) 

m>l,n>l 



where the overall factor —2yG 5(y — L) is omitted, and we have suppressed the or- 
der one coupling constant in front of each operator. The contributions of the above 
operators can be understood easily from symmetry. From the Appendix, in the triplet 
representations, the representation matrix for each element of the Klein four subgroup 
K 4 = {1, Gi = TSTS 2 , G 2 = TST, G 3 = S 2 } is 



G 3 = \ -10, (43) 






where the sign " + " for the representation 3 and " — " for 3'. It is straightforward to 
check that 

Gi(x) = G 2 ( X ) = G 3 ( X ) = (x) , 
-GM = G 2 (cp) = -GM = {<p) . (44) 

As a result, we have 

Gi ((X"V 2 ") 3 > = G 2 ((x"V 2n ) 3 > = G 3 ((x"V 2 ") 3 > = ((X"V 2 ") 3 > , 
-G 1 ( (xV* +1 ) 3( 3 > = G 2 ( (x"V 2n+1 ) 3(30 > = -G 3 ( (x m V 2n+ % {3/) > = < (xV n+1 )3(3 > 

Thus, the vacuum alignments of (x"V 2n ) 3 an d the combination (x"V 2n+1 ) 3 3 ' are 

((xV n ) 3 >oc(0,0,l) T , <(xV n+ %}°c(0,l,0f, <(xV n+1 ) 3 )oc(l,0,0f. (45) 

They are along the same directions as the VEVs (x), (<p) and ((xv 9 ) 3 )5 respectively As 
a result, the high dimensional terms in Eq. ( 1421 only lead to corrections to the electron, 
muon and tau masses, all the off-diagonal entries of the charged lepton mass matrix are still 
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vanishing, and their contributions can be absorbed by a redefinition of the LO parameters 
y e , y^ and y T . This means that the charged lepton mass matrix is diagonal even when 
higher order operators are taken into account, and no rotation of the charged lepton fields 
to the mass eigenstates is needed. 

On the UV brane, the NLO corrections to the Dirac neutrino Lagrangian of Eq. (j2"3"|) 
can be written as 

5C UV = -^Tr (£#6) (p + ^Tr (£#&) ( P 0) 3 + -^Tr (£ x if£ 2 ) (p<f>) s , + h.c. , (46) 

where the overall factor y2G 5(y) is neglected. Expanding the above terms according to 
the contraction rules listed in the Appendix, and after the £4 flavor symmetry breaking 
by the VEVs in Eq. (1201) . we obtain the corrections to the neutrino mass matrix 

/ ei -V2e 3 e 2 -e 3 \ 
5m v = V2 e 3 ei -^ , (47) 

Ve 2 + e 3 -2e ± / V2 

where the parameters £1,2,3 are given by 

6 1 = (A ^3/2 Ya' F ^ = > & °h Cu) » 

e 2 = v 7 ^ ^ ^ /2 ^/ F (y = 0, 0, q, c) , 



(AL) 3/2 A 2 

(AL) 3/2 A 2 



e 3 = V2-^—^ F(y = 0, (3, q, c) , (46 



which are suppressed by v p /A with respect to the LO light neutrino parameters a and b 
of Eq. (|30|) . To first order in the small parameters 6j(z = 1, 2, 3), the lepton mixing angles 
are modified to 



sin 



1 



l.-i 



a/2 a (a& — ofo) ei + a|o| 2 (e 2 — 3e 3 ) 



\/3 [3|6| 4 +(ao-ao) 2 ] 

+V2 [(a 2 + 3b 2 )b - \a\ 2 b] e x + ab 2 (e 2 + 3e 3 

sin 2 0i2 = - H j= — 7-= r--| oei - \/2 6e 2 + c.c. 

2 2A/3 z (ab - ab) I 

S1 " 2 *■ = 1 " apnfc^r ^ (aF + #|2) £I " |6|2Sf2 

- [(2a 2 + 3b 2 ) b - 2\a\ 2 b] e 3 + c.cl , (49) 

where a bar on a letter indicates complex conjugation. It is clear that the NLO corrections 
to the three mixing angles are of order v p /A. We see that, in order to be compatible with 
the current experimental values of the mixing angles [3TH33] . we need a parameter v p /A of 
the order 0.1 ~ 0.2, which is approximately the size of the Cabibbo angle A c . In particular, 
the resulting reactor mixing angle 0i 3 would be expected to be of the same order if no 
accidental cancellations among the coefficients occur. Regarding the light neutrino masses, 
these corrections introduce deviations in e 2 (i = 1, 2, 3) with respect to the LO predictions 

16 



of Eq. (IM|) . and therefore the neutrino spectrum should also be an inverted hierarchy even 
after the high dimensional terms are taken into account. 

In order to see in a quantitative way the behavior of the lepton mixing parameters under 
the NLO corrections, we perform a numerical analysis in the following. The parameter 
V/A has been fixed at the representative value of 0.15. where V stands for the VEV of 
any flavon field in the model. All the undetermined order one dimensionless couplings 
are randomly generated complex numbers, the absolute values are taken to be uniformly 
distributed in the interval [1/2,2], and the phases are chosen to be flatly distributed in 
[0,27r). Then we calculate the associated neutrino masses and the mixing angles which 
are required to lie in the 3a ranges [3TH33] . With the procedure described above, we find 
that the neutrino mass spectrum is inverted hierarchy for all the generated points, which is 
consistent with the conclusion reached from the above analytical estimates. Determination 
of the neutrino mass hierarchy is one of the primary physics goals of the next generation 
precise neutrino oscillation experiments. If the mass spectrum is determined to be normal 
hierarchy in future, our construction would be ruled out. The distributions of the three 
mixing angles and the CP phase 5 are plotted in Fig. |3j From this figure, we can hardly 
see any specifically preferred pattern for 612, and smaller #13 is slightly favored as #13 = 
at LO. The distribution of 5 is essentially random, since our model doesn't introduce 
any mechanisms to predict certain specific values. However, the atmospheric neutrino 
mixing angle #23 with sin 2 #23 > 1/2; which means the second octant of #23? is clearly 
preferred. This fact can be easily understood as follows. The atmospheric mixing angle 
is predicted to be the DC value sin 2 9°^ = 2/3 > 1/2 at LO, after including the NLO 
contributions, the first octant #23 could only be achieved if the NLO corrections are rather 
large. Nevertheless, this requires an additional fine-tuning of the parameters which has 
been reproduced in our numerical simulation only partially and by a few points. At 
present, significant deviation of #23 from maximal mixing value 45° is favored at the level 
of 1.7 — 2 a [3TH33] . and we still don't know which octant #23 lies in. However, future 
long-baseline neutrino oscillation experiments open up the possibility of determining the 
octant of #23 [H]- If future experiments discover that #23 belongs to the second octant 
and the deviation from maximal mixing is somewhat large, the present suggestion of DC 
mixing as a starting mixing pattern of model construction would be strongly favored. 

The deviation of #23 from maximality, if confirmed in future, would have profound 
implications for neutrino mass models based on flavor symmetries. Since the DC mixing 
is preferred for the second octant of #23? we would like to propose another new texture 
which could be taken as a LO approximation for the first octant of #23- By permuting the 
second and third rows of the DC mixing matrix, we find a new interesting mixing pattern 
which is called "modified" DC (MDC) mixing with 



U 



MDC 



/ -L- L. n \ 

111 
v 7 ^ v^ \/3 

, _J L /I 1 

\ V6 v/6 V 3 / 



(50) 



It gives rise to the lepton mixing angles 

sin 2 9% DC = 0, sin 2 9^ DC = 1/2, sin 2 9% DC = 1/3 . (51) 
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Figure 3: The probability distribution of the lepton mixing parameters sin 812, sin #13, sin 823 and the 
Dirac CP phase S in the DC model. 



MDC 
23 



7T 



/2 — #^ c '. Similar to the DC mixing case, the MDC pattern 



Obviously we have 9 

should be corrected by terms of 0(X C ) in order to achieve agreement with the experimental 
data. It is easy to modify the DC model of section [3]to get an MDC model by embedding 
(£i>£i\£i) instead of (£i,£i, £1) into an S4 triplet 3 and by exchanging the assignments 
for £3 and £| under the flavor symmetry. Then we can straightforwardly check that the 
charged lepton mass matrix is also diagonal even in the presence of the high dimensional 
operators. The LO and NLO neutrino mass matrices m' u and Sm' u can be obtained from 
m u and 5m u respectively by switching the second and the third rows 



m,. 



P 2 3«^, Sm' v = PzaSrriu , (52) 

where m„ and 5m u in Eq. fT29|) and Eq. (I47p are the predictions of the above DC model 
for the LO and NLO neutrino mass matrices respectively, and P23 is a permutation matrix 

/100\ 

P23 = 1 . (53) 

\0 1 0/ 

Note that P23 arises due to the exchange of the order of £^ and £f within the triplet 3. As 
a result, the lepton flavor mixing matrix is of the MDC form at LO, the predictions for the 
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Figure 4: The distribution of the atmospheric mixing angle sin 2 023, where the LO mixing is the MDC 
texture. 



neutrino masses are given in Eq. ( )34l) . and thus the neutrino mass spectrum is inverted 
hierarchy as well. Following the method outlined in the DC case, we perform a numerical 
analysis to account for the NLO corrections, the inverted hierarchy neutrino spectrum 
turns out to be still preserved. The predictions for the distribution of the atmospheric 
mixing angle 623 is displayed in Fig. HJ as expected that the first octant of 823 is obviously 
favored. Regarding the mixing angles #13 and # 12 and the Dirac CP phase 8, no specific 
values are preferred within 3a, hence not shown in the figure. 



5 Electroweak constraints 

The electroweak precision measurements (EWPMs) provide stringent constraints to 
new physics at the TeV scale. It is well-known that the custodial symmetry can relax the 
bound of the KK scale [6] in the warped extra dimension. It has been shown that an A 4 
lepton model with an IR brane localized Higgs [2D] can pass the EWPM constraints for 
the KK scale of order 3 TeV if fermions have proper bulk masses. In the following, we 
shall extend the analysis to the bulk Higgs case in the warped extra dimension with the 
custodial symmetry SU(2)l x SU(2)r x U(1)x x Plr, and we shall check that our model 
can indeed pass the EWPM constraints. In addition, we find it is easier for a bulk Higgs 
to pass the EWPM constraints than a brane Higgs. 

We introduce the 4D gauge coupling constants g and g', which are related to the 5D 
gauge coupling constants via g = g$/VL and g' = g' 5 /y/L. A simple way of checking the 
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EWPM constraints is to canonically normalize the SM gauge fields and to determine the 
parameters g, g', and v by requiring that the well-measured values of the Z boson mass 
m z , the fermi constant G F , and the electrical charge e (or the fine-structure constant a) 
be reproduced. Then the corrections to the EW precision observables would be contained 
in the gauge couplings of the fermions. We will use the notation that the SM parameters 
are written with a hat on top of them. The SM expressions for e, Gf and mz in terms of 
the g, g' , and v are given by 



99' 



'g2 + g> 2 



2 I "2 , ~ /2 \ ~ 2 

m z = 4(0 +9 )v 



G, 



1 



V2 



(54) 



In the warped extra dimension, they are corrected by the massive KK modes, and the 
main contributions originate from the gauge sector. Following [57], we define the zero- 
momentum gauge boson propagators for the massive KK modes with BCs (++) and (— +), 
respectively: 



G++(y,y') 

G- + (y,y') 



Ak(kL) 
1 



,2fcL 



,2k y< 



2k 
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where y < (?/>) is the minimum (maximum) of y and y' . To the first order of v 2 /k' 2 , the 
parameters e, Gf, and mz in the present warped model are expressed as [57J 
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where c = g/yg 2 + g' 2 and s = g'/yg 2 + g' 2 - The dimensionless parameters 5± + and 
G l ± + are corrections induced by the exchange of massive KK gauge bosons with BCs 
(±+), and G^ + represents additional KK exchange corrections to 4-fermion interactions 
that contribute to muon decay. They are explicitly given by 
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(57) 



Imposing the match conditions in Eqs. (154)) and 156]) . we obtain the relations between g, 
g', v and g, g', v in the leading order of v 2 /k' 2 
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where oL and d v parameterize the deviations from the SM gauge couplings and Higgs VEV 



nn 



9 - s 2 2k' 2 



[s 2 (5l + -5l + )+c*(2G l ++ -2G l _ + -Gr + )) 



d, 



"1 2 

cm z 
2k' 2 



K + + ^ 2 



2G\ 



2G l 



Gf + ] 



(59) 



respectively. Here s (c) is the sine (cosine) function of the Weinberg angle in the SM 

e 2 / (V2G Fm 2 z ) 2. We can see that the 5D Higgs VEV depends 



with s 2 = 1/2 - J\ 

on the fermion bulk mass q. When the SM SU{2) L doublet is localized close to the UV 
brane, which is required to generate the charged lepton mass hierarchy, the Higgs VEV 
is essentially independent of q. Given the SM Higgs VEV v ~ 246 GeV, we arrive at 
v ~ 250 GeV. Obviously this Higgs VEV is slightly above the SM value 246 GeV. It is 
generally true in the warped extra dimension due to the suppression of the Z gauge boson 
profile near the IR brane and the additional contributions to the gauge boson masses from 
the wavefunction curvature terms [55] . Then we turn to discussing the Z boson couplings 
to the LH charged leptons and RH charged leptons, which are given by 
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for a = e, /i, r, respectively. Since the RH charged leptons have different bulk masses, 
their couplings g^ are flavor- dependent. We introduce parameters d L = (g L — g L )/g L and 
d'x = {.gR — Qp^/gR to measure the deviations of the Z-lepton couplings from the SM values 
g L = g(—l/2 + s 2 )/c and g R = gs 2 /c. Inserting Eq. ( 158]) to Eq. f )60|) . we obtain d L and 
d R in terms of the SM parameters and k' as 
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The charged lepton couplings to the Z boson have been measured precisely to the mille 
level [72]. We require that both deviations \di\ and \d R \ be smaller than 2 x 10 -3 [20] . 
Figure \5\ shows the variation of dp and d R as a function of q and c a , respectively. The 
results of both the bulk Higgs and brane Higgs scenarios are presented. For the bulk 
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Higgs, its mass parameter (3 = 0, 0.25, 0.5 has been taken. We see that the bulk Higgs 

scenario is less stringently constrained than the brane Higgs, and the similar property was 

found in the minimal RS model [9]. In most cases, the deviations are sufficiently small 

and decrease with the increase of /?. Both q and c a are allowed to vary in a large range. 

In particular, all the parameter space for fermions localized near the UV brane (q > 0.5 

and c a < —0.5) can pass the EWPM constraints. 

It is remarkable that the deviation parameters d^ and d% are almost independent of 

the bulk masses of charged leptons, if they are localized close to the UV brane. These 

results can be understood analytically as follows. Since q > 0.5 and c a < —0.5 are valid 

for UV localized charged leptons, ignoring the terms with the small factors e^ 2ci ' kL and 
e (i+2c a )kL^ we bt a }ri 
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Although G l ±+ , G± + and G++ are functions of fermion bulk masses generally, they become 
essentially independent of fermion profiles in this case [57], and Eq. ( 162]) is simplified to 
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(64) 



Obviously d^ and d% have negative and positive contributions, respectively. For the generic 
scale k' = 1.5 TeV in the custodial warped extra dimension, we obtain that both magni- 
tudes are always smaller than the EWPM bounds 2 x 10~ 3 , as shown in Fig. [5j 

Furthermore, the model should also be constrained by the lepton flavor violation pro- 
cesses, in particularly the radiative fi — > er) decay which appears at one-loop level. Due to 
the protection of both the S4 flavor symmetry and the custodial symmetry, we expect that 
our model should more easily satisfy this type of constraints than the conventional warped 
models without flavor protection mechanisms. However, a precise and detailed calculation 
of the lepton flavor violations requires much work and will be presented elsewhere. 

6 Leptonic mixing from residual symmetries within 
S4 flavor symmetry 

In the theory of flavor symmetry, the full Lagrangian is invariant under a family sym- 
metry group Gf at high energies above the symmetry breaking scale. Then Gf is spon- 
taneously or explicitly broken into its subgroups G u and Gi in the neutrino and charged 
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Figure 5: Deviations from the SM Z couplings with the LH charged lepton £l an d the RH charged lepton 
£r in the warped extra dimension with custodial symmetry SU(2)l x SU(2)r x U(1)x X Plr- The solid 
lines are the results for the bulk Higgs scenario considered in the present work, where /? is taken to be 
(green), 0.25 (red) and 0.50 (blue), respectively. The result for the brane Higgs scenario is presented in 
dashed line for comparison. We take a(ft = m z ) = 1/128, m z = 91.1876 GeV, G F = 1.1664xl0~ 5 GeV~ 2 , 
and k' = 1.5 TeV as input parameters, and for d^, we set q = 0.51. The dotted lines stand for the current 
experimental bounds (±2 x 10 -3 ). 



lepton sectors at LO, respectively, and certain lepton mixing matrix Upmns arises as a 
result of the misalignment between the two residual symmetry groups G u and Gi [73] . In 
general, G u and G\ should be different subgroups of Gf, otherwise there would be no mix- 
ing among the three mass eigenstates. We embed the LH lepton doublets into a unitary 
three-dimensional representation p of Gf to generate a realistic lepton mixing pattern [71] . 
The invariance under the residual symmetries G u and Gi requires 

p 1 '(g v )m u mlp(g v ) = m v m\, p ] (gi)mim\p(gi) = mim], g v G G u , g x 6 G x , (65) 

where the neutrino mass matrix m v and charged lepton mass matrix mi are defined by 



C 



-l L mil 



r ~ VLm u v R 



(66) 
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'rm- Here neutrinos 



and l L = (e, p, r)£, l R = (e, p, r)%, v L = (u e , v^ v T )\ and u R = (u e , u t 
are assumed to be Dirac particles. If one chooses G v or Gi to be a non-Abelian subgroup, 
according to the Schur's theorem, this would result in a complete or partial degeneracy 
of the mass spectrum, which is not compatible with the observation that both the light 
neutrinos and the charged leptons have different masses. Therefore we shall restrict G v 
and Gi to be Abelian subgroups. The representation matrices p(g u ) and p(g e ) can be 
diagonalized by two unitary transformations U v and U e as follows 



Uip{g v )U u = p\g u ), U}p{gi)Ui = p(gi 



(67) 



where p(g u ) and p(gi) are diagonal unitary matrices, and g v and gi can be taken to be the 
generators of the Abelian subgroups of G v and Gi, respectively. We note that the matrices 
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U u and U[ are determined uniquely up to diagonal unitary matrices K v j and permutation 
matrices P v j, respectively 

U v ,i -> U^K^P^ . (68) 

From the invariance requirement of Eq. (165]) , it follows 

Ulm v m\U v = U ] v p ] (g p )m p mlp(g u )U p = p (g v )Ulm u mlU v p(g u ) , 
U{mim\Ui = U{ p ] \gi)mim\p(gi)Ui = p*(gi)U}mim}Uip(gi) . (69) 

As a result, both Ulm v m\U v and [//m^mJC// have to be diagonal, i.e., the hermitian combi- 
nation m v m\ is diagonalized by U u , and mim] is diagonalized by Ui as well. Consequently, 
we can determine the lepton mixing matrix Upmns as 

Upmns = U}U V (70) 

up to permutations of rows and columns. Obviously if G\ and G v are the same subgroups of 
5*4, the PMNS matrix would be an identity matrix, which is incompatible with the present 
data [3TH33] . In the following, we perform a comprehensive study of the possible lepton 
mixing patterns within S4 flavor symmetry by scanning the different choices for G v and 
G[. The non-trivial Abelian subgroups of S4 are Z 2 , Z 3 , Z4 and K$, where K A = Z 2 x Z 2 
is the Klein four group. The Z 2 subgroups are irrelevant to us, since at least two of the 
eigenvalues are degenerate, such that the three generations of neutrinos or charged leptons 
can not be distinguished. On the technique side, this degeneracy prevents us from uniquely 
pinning down the lepton mixing matrix, because the unitary transformations U v and C/j 
can not be fixed unambiguously in this case. The group structure of S4 has been discussed 
in Ref. [75] in detail, and concretely it has four Z 3 subgroups, three Z 4 subgroups and four 
K 4 subgroups. In terms of the generators S and T, they can be expressed as 



Z 3 subgroups 



ZW = {1,T,T 2 }, zf = {l,T 2 5 2 ,5 2 T}, 

Zf = {1, S 2 TS 2 , STS} , Zf = {1, S 2 T 2 , TS 2 } . 



(71) 



• Z4 subgroups 

Z[ 1] = {l,S,S 2 ,S 3 }, Zf = {l,TS,T 2 S 2 T,T 2 ST}, zf = {l,ST,TS 2 T 2 ,TST 2 } 



K4 subgroups 



K[ l) = {1, STS 2 , T 2 S, TS 2 T 2 } , Kf ) = {1, TSTS 2 , TST, S 2 } , 

Kf ] = {1, ST 2 , S 2 TS, T 2 S 2 T} , Kf = {1, TS 2 T 2 , S 2 , T 2 S 2 T} . 



(72) 



In the following, all the possible Abelian remnant symmetry groups G v and G\ of £4 are 
considered one by one, the corresponding lepton mixing matrix Upmns and the group 
structure generated by G v and G\ are listed in Table |2j Since we are mainly interested 
in the mixing angles, we only present the absolute values of the mixing matrix entries, 
which are denoted by ||C/pmats||- From this table, we see that we have a broader set of 
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remnant symmetry groups Gi and G v for Dirac neutrinos, recalling that G v is generally 
fixed to K4 in the case of Majorana neutrinos. Clearly seven leptonic mixing patterns 
up to permutations of rows and columns can be produced within 5*4, and only four of 
them remain if we require the elements of Gi and G v generate the whole group £4. It is 
interesting to note that if the S4 flavor symmetry is broken into G v = Z4, Gi = Z 3 or 
G v = Z 3 , Gi = Z&, where Z3 and Z± denote any Z3 and Z4 subgroups, respectively, we will 
be able to produce the so-called TFH mixing pattern [31H36] . The same point has also 
been observed in Ref. [M] |j. However, if neutrinos are Majorana particles, to generate 
TFH mixing without fine tuning, the flavor symmetry group should be A (96) or any group 
containing it. This implies that, in contrast to A(600), {Z\% x Zq) x S3 and A(1536) for 
Majorana neutrinos, smaller family symmetry groups could possibly yield lepton mixing 
patterns which lie within the 3a ranges of the current global fitting values, if neutrinos 
are Dirac particles [76]. Similar to the Majorana neutrinos case, the combinations of 
G v = K\ ' ' and G\ = Z 3 lead to the tri-bimaximal mixing. On the other hand, we are 
able to derive the DC mixing for G u = Z 3 and Gi = K\ ' , which is exactly the LO 
symmetry breaking chain of our model in section [31 Finally, the bimaximal mixing can be 
produced if the residual symmetry groups in the neutrino and charged lepton sectors are 
K4 or Z4, please refer to Table [2] for detailed possible choices of G v and G\. 



7 Summary and conclusions 

The remarkable discovery of a rather large mixing angle #13 and the recent indication 
of a sizable deviation of the atmospheric neutrino mixing angle #23 from 7r/4 have pro- 
found impact on neutrino physics. Many neutrino mass models are ruled out due to their 
predictions of small #13. In the present work, from the viewpoint of flavor symmetry, we 
have investigated the proposal that the LO neutrino mixing matrix is the DC pattern and 
then it is corrected by the NLO contributions of order A c such that the resulting lepton 
mixing angles can be in the experimentally preferred ranges. 

It is well-known that warped extra dimensions provide an attractive solution to the 
gauge hierarchy problem, an understanding of the observed hierarchies in the fermion 
masses and mixing angles, and a natural suppression of the dangerous flavor changing 
neutral current processes. Inspired by the successes of discrete flavor symmetry in pre- 
dicting the lepton mixing angles, we combine the warped extra dimension and the discrete 
flavor symmetry together in this work. We construct a lepton model based on the fla- 
vor symmetry S4 x Z2 x Z<i in the warped extra dimension with the custodial symmetry 
SU(2)l x SU(2) R x U(l)x x P LR . In this model, all the matter fields, gauge fields and 
the Higgs field are allowed to propagate in the bulk. The flavon fields breaking the S4 
flavor symmetry in the neutrino and charged lepton sectors live on the UV and IR branes, 
respectively. The light neutrinos are taken to be Dirac particles, the correct neutrino mass 
scale can be obtained by choosing the fermion and Higgs bulk mass parameters, and the 



9 Only the choice G v = Z4 and G; = Z 3 was pointed out in Ref. 
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Table 2: The leptonic mixing patterns and the generated group structures for the S4 flavor symmetry 
breaking into different subgroups G v and Gi in the neutrino and charged lepton sectors, respectively, 
where Z3 (Z4) denotes any Z3 (Z4) subgroup. D4 stands for the dihedral group of degree 4 and order 8, 
and geometrically it is the symmetry group of the square. 



charged lepton mass hierarchy is generated as usual in RS models via wavefunction over- 
laps. The lepton mixing matrix turns out to be of the DC form at LO, and the neutrino 
mass spectrum is predicted to be inverted hierarchy. 

For DC mixing under the hypothesis of Dirac neutrinos, although the symmetry group 
of the left-handed neutrino fields is Q v = U(l) x 17(1) x 17(1) x U(l) in the charged 
lepton diagonal basis, it is sufficient to produce DC mixing by just requiring that the 
neutrino sector be invariant under the action of one element h E Q v with non-degenerate 
eigenvalues o Therefore the DC mixing can be generated if the flavor symmetry group 
Gf contains the element h which can be preserved in the neutrino sector. This is the 
reason why DC mixing can be derived from the S4 flavor symmetry, and h corresponds to 
the generator T in the present work, as shown in section [2j 



id 



This statement is generally applicable to the mass-independent textures in the case of Dirac neutrinos. 
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Although the exact DC mixing leads to the mixing angles sin 2 8y 3 c = 0, sin 2 8® 2 C = 1/2 
and sin 2 9%f = 2/3, which is excluded by the present data, the LO predictions are modified 
by the high dimensional operators allowed by the symmetries of the model. The induced 
deviations of all the three mixing angles from the DC values are linear in the expansion 
parameter V/A, where V denotes the VEV of any flavon field. It is usually assumed that all 
the VEVs are approximately of the same order of magnitude. In order to achieve agreement 
with the experimental measurements, the parameter V/A should be of order A c . Detailed 
numerical simulations show that the neutrino mass spectrum is still inverted hierarchy and 
the second octant of 823 is preferred after the NLO contributions are taken into account. 
Therefore our model can be tested in next generation neutrino oscillation experiments, 
which are deliberately designed to resolve the neutrino mass ordering, the octant of 823 
and eventually determine the leptonic CP violating phases. Since the first octant of 823 is 
also allowed by the present data, we suggest that the "modified" DC (MDC) mixing as a 
LO approximation and subsequent 0(\ c ) corrections can account for the data on neutrino 
mixing for the case of 823 in the first octant, where the MDC mixing matrix can be obtained 
by permuting the second and the third rows of the DC mixing matrix. By exchanging the 
assignments for £| and £f and unifying (£i)6i>£i) instead of (£i,£i, £1) into a £4 triplet 
3, our constructed DC model would become a model for the MDC mixing. Moreover, we 
study the constraints on the model from the electroweak precision measurements, which 
are estimated by requiring small deviations from the SM tree-level couplings. We show 
that the model is compatible with the experimental constraints for usual KK mass scales 
of order 3 TeV. 

Finally we have investigated the associated lepton mixing patterns for all the possible 
remnant subgroups G v and G\ in the neutrino and charged lepton sectors respectively 
within S4, where neutrinos are assumed to be Dirac particles. It is notable that the TFH 
mixing can be generated within £4 instead of A (96) for Dirac neutrinos if the S4 flavor 
symmetry is broken into G u = Z±, Gi = Z 3 or G v = Z 3 , Gi = Z4. We can derive the DC 
mixing from the combinations of G u = Z3 and Gi = JQ ' . This is the guiding principle 
of our model building. In addition, the tri-bimaximal mixing and bimaximal mixing can 
also be produced if we require the elements of G v and G\ give rise to the entire group £4. 
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Table 3: Character table of S4 group, where hd denotes the order of the elements of Cj, and rid is the 
number of the elements in the class C,- . 



Appendix: The group theory of S± 

Si is a symmetric group of degree four. Geometrically, it is the group of orientation- 
preserving symmetries of the cube (or equivalently, the octahedron). The group has 24 
distinct elements, and it can be generated by two generators S and T which satisfy the 
relations 



,s 4 = r 3 = i, ST 2 S = T. 



Without loss of generality, we could choose 

S = (1234), 



(123) 



(73) 



(74) 



where the cycle (1234) denotes the permutation (1,2,3,4) — > (2,3,4,1), and (123) means 
(1,2,3,4) — > (2,3,1,4). The 24 group elements are divided into 5 conjugacy classes as 
follows: 



C 2 

c 3 
c 4 



STS 2 = (12), TSTS 2 = (13), ST 2 = (14), S 2 TS = (23), TST = (24), T 2 S = (34) 

TS 2 T 2 = (12) (34), S 2 = (13)(24), T 2 S 2 T = (14)(23) 

T = (123), T 2 = (132), T 2 S 2 = (124), S 2 T = (142), S 2 TS 2 = (134), STS = (143), 

S 2 T 2 = (234), TS 2 = (243) 

S = (1234), T 2 ST = (1243), ST = (1324), TS = (1342), TST 2 = (1423), 

S 3 = (1432) 



Since the number of the irreducible representations is equal to the number of conjugacy 
classes, S4 has five irreducible representations: two singlets 1 and 1', one doublet 2, and 
two triplets 3 and 3'. We note that both 3 and 3' are faithful representations. The 
character table of S4 is presented in Table El and the Kronecker products between two 
irreducible representations, which can be easily deduced from the character table, are 
listed below: 

1®R = R®1 = R, 1'®1' = 1, l'®2 = 2, l'<8>3 = 3', 1' <8> 3' = 3, 

2 <g> 2 = 1 © 1' ©2, 2©3 = 3©3', 2 © 3' = 3 © 3', 3©3 = 1©2©3©3', 

3 © 3' = 1' © 2 © 3 © 3', 3' © 3' = 1 © 2 © 3 © 3' , (75) 
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Table 4: Representation matrices for the generators S and T in different S4 irreducible representations. 



where R denotes any £4 irreducible representation. Different presentations of the S4 group 
have been discussed in literatures [75 |[77H80] . In the present paper, we would like to work 
in the charged lepton diagonal basis such that the DC mixing completely comes from the 
neutrino sector, the representation matrices of the generators S and T for different S4 
irreducible representations are presented in Table HI The adopted representation here is 
related to other existing choices by unitary transformations. In the following, we list the 
Clebsch-Gordan coefficients derived from this basis. All the Clebsch-Gordan coefficients 
are reported in the form ofa<8>/3, where the «; are the base vectors of the representation on 
the left of the product, and the f3j are from the representation on the right of the product. 

• 1' ® 2 = 2 



1' <8> 3 = 3' 



1' (8) 3' = 3 



3' 



-ai/?i 




2(8)2 



l'©2 




1 ~ «i/3i + a 2 (3 2 
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• 2 © 3 = 3 © 3' 



• 2 © 3' = 3 © 3' 



3©3 = 1©2©3©3' 



1' ~ «i/3 2 - a 2 /3± 
a 2 /3 2 - ai/3i 



2 

«i/3 2 + a 2 /3i 



ai/3i + A/3a 2 /3 2 

ol\$i + A/3a 2 /3i 

-2ai/3 3 

\/3ai/3i - a 2 /3 2 
2a 2 (3 3 



y/3ai(3 2 - a 2 fii 
3 ~ \/3ai/?i - a 2 /3 2 
2« 2 /3 3 

aiA + •v / 3a 2/ 5 2 
3' r | «i/?2 + v / 3a 2 /3 i 
-2ai/3 3 



1 ~ ai/3i + a 2 /3 2 + a 3 /3 3 



^ , «i/3i + a 2 /3 2 - 2a 3 /3 3 
v^ai/32 + a/3q; 2 /3i 



«l/?3 + «3/?l 

[ -a 2 (3 3 - a 3 (3 2 
«i/3i - a 2 (3 2 

' a 2 (3 3 - a 3 /3 2 

| a 3 /?i - «i/3 3 

ai/3 2 -a 2 /3i 
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• 3 © 3' = 1' © 2 © 3 © 3' 



• 3' (8) 3' = 1 © 2 © 3 © 3' 



1' ~ «i/3i + a 2 /3 2 + a 3 (3 3 
-V3a l (3 2 - V3a 2 fti 



2 

«i/3i + a 2 /3 2 - 2a 3 /3 3 



' «2/?3 - "3/32 

| a 3 /3i - ax/33 
a t f3 2 - a 2 Pi 

anfo + aizPi 

| -a 2 (3 3 - a 3 (3 2 

ai/3x - a 2 (5 2 



1 ~ ax/?! + a 2 /3 2 + a 3/ g 3 



^ , «i/3i + a 2 /3 2 - 2a 3 /3 3 



ai/?3 + a 3 /3i 
[ -ct 2 /3 3 - a 3 /3 2 
ai/3i - a 2 /3 2 

' a 2 (3 3 - a 3 (3 2 

| a 3 /3i - «i/3 3 

.«l/?2 -«2/3i 
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